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Approximate solutions for the gluon and ghost propagators as well as the running coupling 
in Landau gauge Yang-Mills theories are presented. These propagators obtained from the 
corresponding Dyson-Schwinger equations are in remarkable agreement with those of re- 
cent lattice calculations. The resulting running coupling possesses an infrared fixed point, 
Q s(0) = 8.92/A r c for all gauge groups SU(JV C ). Above one GeV the running coupling 
rapidly approaches its perturbative form. 



The investigation of Quantum Chromo Dynamics (QCD) over the last decades made obvious that 
a description of hadrons and their processes from the dynamics of confined quarks and gluons 
in one coherent approach is an outstanding complicated task. Knowledge of a completely non- 
perturbative running coupling would be a very important step to make towards such a QCD-based 
description of hadrons. 

Given the experimental results in hadron physics it is evident that baryons and mesons are 
not elementary particles in the naive sense of the word "elementary". The partonic substructure 
of the nucleon has been determined to an enormous precision leaving no doubt that the parton 
picture emerges from quarks and gluons, the elementary fields of QCD. This contrasts the well- 
known fact that these quarks and gluons have not been detected outside hadrons. This puzzle 
was given a name: confinement. The confinement hypothesis was formulated several decades 
ago, nevertheless, our understanding of the confinement mechanism(s) is still not satisfactory. 
Moreover, in contrast to other non-perturbative phenomena in QCD {e.g., dynamical breaking of 
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1 Motivation: Some aspects of confinement 
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chiral symmetry, Ua(1) anomaly, and formation of relativistic bound states), it seems not even 
clear, at present, whether the phenomenon of confinement is at all compatible with a description 
of quark and gluon correlations in terms of local fields in the usual sense of quantum field theory. 

One important feature of QCD is asymptotic freedom: Employing the renormalisation group 
(RG) in perturbative calculations clearly leads to an understanding of the experimental observa- 
tion that partons interact at large (spacelike) momentum transfer only weakly. On the other hand, 
RG based arguments prove that perturbation theory is insufficient to account for confinement in 
four-dimensional field theories: Confinement requires the dynamical generation of a physical 
mass scale. In presence of such a mass scale, however, the RG equations imply the existence of 
essential singularities in physical quantities, such as the S'-matrix, as functions of the coupling 
at g = 0. This is because the dependence of the RG invariant confinement scale on the coupling 
and the renormalisation scale fi near the ultraviolet fixed point is determined by [1] 

A = MeXP (-/^))" %eXP (-^)' /30> °- (1) 

Therefore a study of the infrared behaviour of QCD amplitudes requires non-perturbative meth- 
ods. In addition, as infrared singularities are anticipated, a formulation in the continuum is 
desirable. One promising approach to non-perturbative phenomena in QCD is provided by stud- 
ies of truncated systems of its Dyson-Schwinger equations, the equations of motion for QCD 
Green's functions, for recent reviews see e.g. [2]. As we will see in the following these studies 
also allow to extract a non-perturbative running coupling. One word of warning is, however, in 
order: A non-perturbative running coupling is not a uniquely defined object. The extension of the 
running coupling into the infrared domain requires in the first step its definition from a specifi- 
cally chosen Green's function. Of course, within a chosen gauge the result should be unique, and 
one should be able to prove this uniqueness from the Slavnov-Taylor identites of QCD. Noting 
that Green's function are not gauge invariant it is not at all obvious whether the comparison of 
non-perturbative running couplings defined in different gauges can be meaningful at all. 



2 A possible definition of the non-perturbative running coupling in Landau gauge 

As stated above, the definition of the non-perturbative running coupling rests on a specifically 
chosen Green's function. To this end we note that the ghost-gluon vertex in Landau gauge ac- 
quires no independent renormalisation, in specialist language Z\ = 1. This relates the charge 
renormalisation constant Z g to the ones for the gluon and ghost wave functions, 1 = Z\ = 
Zgy/Z^Zs: The gluon leg provides a factor \fZ~3, the two ghost legs Z3. As we will demon- 
strate in the following this allows for a definition of the running coupling resting solely on the 
properties of the gluon and ghost propagators. 

In linear covariant gauges the gluon propagator is of the form 5 

Z(k 2 ) / k^k u \ k^K 

U^u Gluon - — p— I <V - -p- I +4-^- U) 

The Lorentz condition d^A^ = is strictly implemented only in the limit £ — > which defines 
the Landau gauge. Then the gluon propagator is strictly transverse with respect to the gluon 

5 As usual in these studies a Wick rotation to Euclidean space has been employed. 
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momentum k. The gluon dressing, also in the non-perturbative domain, is solely described by 
the function Z(k 2 ). The Faddeev-Popov ghosts are introduced in the process of quantization: 
The functional integral over these ghosts is a representation of the Jacobian factor induced in the 
generating functional when enforcing the gauge condition. The scalar ghost fields belong to the 
trivial representation of the connected part of the Lorentz group. As local fields with space-like 
anti-commutativity, they violate the spin-statistics theorem and are thus necessarily unphysical. 
The general form of their propagator reads 

£>Ghost = — — fc2~~^' @) 

A comparison to the tree level form of these propagators reveals that Z(k 2 ) — > const, and 
G(k 2 ) — > const., up to the perturbative logarithms, for asymptotically large momenta k 2 . 

In the next step we discuss the employed non-perturbative subtraction scheme and its relation 
to the definition of the running coupling. As already stated, the starting point is the following 
identity for the renormalisation constants 

Z x = Z g Zl /2 Zz = 1 , (4) 

which holds in Landau gauge. It follows that the product g 2 Z (k 2 )G 2 (k 2 ) is RG invariant. In 
absence of any dimensionful parameter this (dimensionless) product is therefore a function of 
the running coupling g, 

g 2 Z(k 2 )G 2 (k 2 ) = f(g 2 (t k ,g)) , t k = \\nk 2 /n 2 . (5) 

Here, the running coupling g(t, g) is the solution of d/dt g(t, g) — (3(g) with g(0, g) = g and the 
Callan-Symanzik /3-function (3(g) = —(3og 3 + 0(g 5 ). The perturbative momentum subtraction 
scheme is asymptotically defined by f(x) — > x for x — > 0. This is realized by independently 
setting 

Z(fi 2 ) = 1 and G(fi 2 ) = 1 (6) 

for some asymptotically large subtraction point k 2 = u 2 . If the quantity g 2 Z (k 2 )G 2 (k 2 ) is to 
have a physical meaning, e.g., in terms of a potential between static colour sources, it should 
be independent under changes (<?, /u) — > (g',fi r ) according to the RG for arbitrary scales fi'. 
Therefore, 

g 2 Z(u> 2 )G 2 (n' 2 ) = g' 2 = g 2 W/n),g) , (7) 

and, f(x) = x, Vx. This can thus be adopted as a physically sensible definition of a non- 
perturbative running coupling in the Landau gauge. In the scheme summarized in the next sec- 
tion, it is not possible to realize f(x) = x by simply extending the perturbative subtraction 
scheme (6) to arbitrary values of the scale /i, as this would imply a relation between the func- 
tions Z and G which is inconsistent with the leading infrared behaviour of the solutions as will 
become evident from the discussion presented in the next section. For two independent functions 
the condition (6) is in general too restrictive to be used for arbitrary subtraction points. Rather, 
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in extending the perturbative subtraction scheme, one is allowed to introduce functions of the 
coupling such that 

Z{iJ?) = f A (g) and G( M 2 ) = f G (g) with f*f A = 1 , (8) 

and the limits j a. g — ► 1 , 9^0. Using this it is straightforward to see that for k 2 ^ [i 2 one 
has(t fe = (lnfc 2 / M 2 )/2), 

Z(k 2 ) = c X p{-2j° {tk ' 9 \l^}f A (g(t k ,g)), (9) 

g(tk,g) 



G(k 2 ) = cxp {-2j 9k ' 9 dl^}f G (g(t k ,g)). 



Here 7.4(3) and 7g(s0 are the anomalous dimensions of gluons and ghosts, respectively, and 
/3(g) is the Callan-Symanzik /3-function. Eq. (4) corresponds to the following identity for these 
scaling functions in Landau gauge: 

27g(<7) + 1a{9) = ~0{g) ■ (10) 

One thus verifies that the product g 2 ZG 2 indeed gives the running coupling (i.e., Eq. (5) with 
f(x) = x). Perturbatively, at one-loop level Eq. (10) is realized separately, i.e., 7g(<?) = 
-60(g)/ 9 and 7.4(3) = -(1 - 25) 0(g) /g with 5 = 9/44 for N f = and arbitrary N c . 
Non-perturbatively one can still separate these contributions from the anomalous dimensions by 
introducing an unknown function e(g), 

lG (g) =: -(6 + e(g))^ j A (g) = -(1 - 25 - 2e(g)) ® . (11) 
9 9 

This allows to rewrite Eqs. (9) as follows: 



z{k 2 } = [«_5^1££j exp|-4y dl-Y^f A (g(t k ,g)), (12) 



g 2 (tk,g) 



s ( rg{tk,g) 



G(k 2 ) = [^f 11 ) c ^{ 2 J g dl^jf G (g(t k ,g)). 

This is also possible in the presence of quarks. In this case one has 5 = / 0o — 9iV c / (447V C — 
8Nf) for Nf flavours in Landau gauge. The above representation of the renormalisation func- 
tions expresses clearly that regardless of possible contributions from the unknown function e(g), 
the resulting exponentials cancel in the product G 2 Z. For a parameterisation of the renormali- 
sation functions, these exponentials can of course be absorbed by a redefinition of the functions 
Ja.g- The only effect of such a redefinition is that the originally scale independent functions 
fA,G(g(tk,g)) will acquire a scale dependence by this, if e 7^ 0. 

For the truncation scheme presented in the next section it is possible, however, to obtain 
explicitly scale independent equations thus showing that the solutions for the renormalisation 
functions G and Z obey one-loop scaling at all scales [3]. In particular, this implies that the 
products g 2S G and g 2 ^ 2 ^ Z are separately RG invariants (as they are at one-loop level). As 
for the renormalisation scale dependence, the non-perturbative nature of the result is therefore 
buried entirely in the result for the running coupling. 
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Fig. 1. Diagrammatic representation of the truncated gluon and ghost Dyson-Schwinger equations studied 
in this letter. Terms with four-gluon vertices have been dismissed. 



Having at hand a non-perturbative definition of the running coupling which employs only the 
properties of the gluon and ghost propagators, an equally non-perturbative method is then re- 
quired to determine these propagators. Their infrared singularity structures are particularly well 
accessible by non-perturbative continuum methods. To this end the Landau gauge Dyson- 
Schwinger equations (DSEs) have been solved analytically in the infrared [3, 4, 5, 6, 7, 8]. 
Necessarily, however, DSE studies will always be subject to truncations in order to obtain the 
closed system of equations to be studied. Thus, the quality of such a study crucially depends on 
the justification of the truncations. Some uncertainty about these truncations remains, however, 
even in the most ambitious study. On the other hand, corresponding lattice calculations, see e.g. 
[9, 10, 11], include all non-perturbative physics but are limited for small momenta by the finite 
lattice volume. Despite their respective shortcomings both these approaches agree encouragingly 
well in the general observations: there is clear evidence for an infrared finite or even vanishing 
gluon propagator and a diverging ghost propagator. This is in accordance with the Kugo-Ojima 
confinement criterion, which in Landau gauge includes the statement that the ghost propagator 
should be more singular than a simple pole [12]. 

For the purpose of this talk we will concentrate on the truncation scheme for DSEs presented 
in refs. [7, 8]. The infrared behaviour of the ghost and gluon system was studied analytically 
for a wide class of non-perturbatively dressed ghost-gluon vertex functions in ref. [6]. Based on 
few and general assumptions about the generic structure of this vertex, it was thereby concluded 
that its dressing should not affect the qualitative findings. We therefore restrict to bare three- 
point vertices for simplicity. All contributions from explicit four-gluon vertices are neglected in 
addition. A diagrammatical representation of the resulting system of equations is presented in 
Fig. 1 . This truncation scheme provides the correct anomalous dimensions of the ghost and gluon 
dressing functions, Z(k 2 ) and G(k 2 ), in the ultraviolet region of momentum. On the other hand, 
it reproduces the infrared exponents found in [5, 6] which are close to the ones extacted from 
lattice calculations [9, 10, 11]. The numerical solution to the truncated DSEs for the gluon and 
the ghost propagators 6 proved to be compatible with only one out of the two solutions reported 
in the infrared analysis of ref. [5]. This thus demonstrates that not every analytical solution for 
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6 A detailed description of corresponding numerical techniques can be found in refs. [13, 4, 7]. 
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DSE vs. lattice results (16 x32) 



DSEvs. lattice results (16 x32) 
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Fig. 2. Solutions of the Dyson-Schwinger equations compared to recent lattice results for two colours [11]. 



asymptotically small momenta necessarily connects to a numerical solution for finite momenta. 

The corresponding infrared behaviour of the propagators is given by I?Giuon(fc 2 ) ~ (/c 2 ) 2k_1 
and L>Ghost(fc ) ~ (fc 2 )^ 1 with n = (93 - \/1201)/98 w 0.595: One obtains a weakly in- 
frared vanishing gluon propagator and a strongly infrared enhanced ghost propagator. In Fig. 2 a 
comparison of the numerical solution of DSEs with the results of recent lattice calculations [11] 
for the case of two colours is shown. As the solutions on the lattice include all non-perturbative 
effects, the results shown in Fig. 2 suggest that the omission of the two-loop diagrams in the 
truncated DSEs mostly effects the region around the bending point at 1 GeV. Given the limita- 
tions of both methods the qualitative and partly even quantitative agreement is remarkable. The 
combined evidence of the two methods points strongly towards an infrared vanishing or finite 
gluon propagator and an infrared singular ghost propagator in Landau gauge. 



4 Infrared fixed point and as (/x 2 ) 

The non-perturbative definition of the running coupling given in section 2 can be summarized as 
follows: 



(13) 



Here we have made explicit the dependence of the propagator functions on the renormalisation 
point. An important point to notice in the results described in the last section is the unique 
relation between the gluon and ghost infrared behaviour. This is no accident: Consistency of the 
DSEs require that the product Z(k 2 )G 2 (k 2 ) goes to a constant in the infrared. 
Correspondingly we find an infrared fixed point of the running coupling: 



a s (0) 



2ir r(3 - 2 K )r(3 + «-,)r(i + k) 



3Nr. 



T 2 (2- k)T(2k) 



93- VT20T 
98 



(14) 



For the gauge group SU(3) the corresponding numerical value is «s(0) ~ 2.972. Of course, this 
result depends on the employed truncation scheme. In ref. [6], assuming the infrared dominance 
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Fig. 3. The strong running coupling (left panel) and the corresponding f3 function (right panel, fully drawn 
line). The latter is compared to the one-loop j3 function (dotted line) and a polynomial in a (dashed line). 



of ghosts it has been shown that the tree-level vertex result as(0) ~ 2.972, among the general 
class of dressed ghost-gluon vertices considered in the infrared, provides the maximal value for 
as(0). If the exponent k is chosen in an interval between 0.5 and 0.7 (as strongly suggested by 
lattice results) one obtains as(0) > 2.5 [6]. 

In the case of three colours values for physical scales can be obtained by requiring the ex- 
perimental value as(Mz = (91.2GeV) 2 ) = 0.118. Together with the numerical solutions for 
the gluon and the ghost propagators we can summarize our knowledge of the running strong 
coupling in the following fit: 

as{x) = Qs(0) — — (15) 

a s (0) = 2.972, o a = 5.292GeV~ 2a2 , a 2 = 2.324, b x = 0.034GeV~ 2b2 , b 2 = 3.169. 

In fig. 3 the running coupling and the corresponding /3 function are shown. 

Finally, we note that indications for an infrared finite coupling have recently also been ob- 
tained within the background field method from the Exact Renormalisation Group Equations 
[14]. 



5 Outlook: a s (p 2 ) and physical observables 



To summarize: We have provided a non-perturbative definition for the running strong coupling 
in the Landau gauge. The underlying picture is related to confinement of transverse gluons, see 
e.g. [15] and references therein. Loosely speaking, one can summarize this by stating that gluons 
are confined by the Faddeev-Popov ghosts. The most interesting result for the running coupling 
is the existence of an infrared fixed point. The occurence of this fixed point can hereby be traced 
back to very general properties of the ghost Dyson-Schwinger equation [6]. Taking only one- 
loop terms in the gluon Dyson-Schwinger equation into account one obtains that cts(/j 2 ) oc 
1/N C with N c being the number of colours. For three colours the corresponding critical value 



8 



R. Alkofer et al. 



for cts{0) slightly depends on the approximations used. Together with the results of lattice 
calculations we confidently conclude that ces(0) ~ 3 or slightly lower. 

Of course, a phenomenological verification of the presented picture would be highly wel- 
come. From so-called analytic pertubation theory, see e.g. [16] and references therein, it is 
known that an infrared finite coupling allows for an almost straightforward calculation of am- 
plitudes relevant in inclusive r decay, in e+-e~-annihilation into hadrons, in inelastic lepton- 
hadron scattering etc. The interesting point hereby is that related observables might give bounds 
on allowed values for as(0). 
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